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D "You haven't measured it, you've only calculated it..." 

Jan Smtt. 

cn 

CN We review nonperturbative calculations of the rate of sphaleron transitions on the lattice in 1 + 1- 

dimensional field theories and introduce a way to perform gauge-invariant Gibbs averages in classical 
non-Abelian Higgs theories. 
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Effective potential of the Chern-Simons variable N in the Standard Model (SM) is 
bounded at infinity. For this reason in hot electroweak plasma one expects Brownian mo- 
CK tion of the system between topologically distinct vacua labelled by integer valines of N : 

A(i) = < (A(i) - A(0))^ > ^ Tt, ast^oo, (1) 

P-i 

^ where ^ • ^ stands for the Gibbs average. Diffusion rate T is the rate of sphaleron transitions. 

The parameter A(t) at large times receives contributions only from the relevant fiuctuations 
interpolating between different classical vacua, which makes it possible in practice to calculate 
r numerically on the lattice. P is a nonperturbative definition of the rate of anomalous 
baryon number nonconservation at high temperatures in the SM [1], which can be used in the 
semiclassical temperature domain and beyond it. At high temperatures one neglects quantum 
corrections and solves classical equations of motion to generate static solutions like kinks [2] 
or sphalerons [3]. 

The correlator (1) was first evaluated in 1 -|- 1-dimensional Abelian Higgs model in [4] and 
the Brownian behaviour of the Chern-Simons number was found. Two methods to evalute 
A(t) are known so far. The first one is based on the microscopical classical real-time Langevin 
equation, derived in [4] in such a way that the white noise would respect gauge invariance. 
The Lagrangian of the model is: 



with = d — iA^. In this Abelian gauge theory one can resolve constraints in the Coulomb 
{diAi = 0) gauge using polar coordinates (j) = /oe'". The only (gauge-invariant!) degree of 
freedom left of the gauge field after the constraints have been resolved is the Chern-Simons 
variable N = l/27r / dxAi. The momentum p„ canonically conjugated to the gauge-depenent 
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degree of freedom a is the the density of electric charge. The white noise cannot be intro- 
duced in the equation for p„, because it would violate the conservation of the electric charge. 
Therefore the noise term must be introduced only to the gauge-invariant variables like p. The 
solution N(t) of the resulting first-order Lagevin equation is obtained numerically and the 
diffusion parameter (1) is evaluated as 

A{t) = — r° dr {N{t+T) - N{T) f , as To ^oo, (3) 

To Jo 

In this way the diffusion (I) was observed [4] with the expected thermal activation behaviour 
of the diffusion rate T ~ exp( — Esp/T) at T < E^p, where E^p — I is the sphaleron mass. 
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Figure I: ln(r) vs (3 for lattice spacings a = (0.16,0.32) in the Abelian Higgs model [5]. 

The second method is to calculate the Gibbs average (I) directly according to the defi- 
nition of the classical Green function: 

A(t) = J VpV(t> exp[-H(p,(t>)/T] Ap^(t) (4) 

where Apif,(t) is evaluated on the microcanonical trajectories p(t), 4'{t) - solutions of the Hamil- 
tonian equations of motion with the initial conditions p{t) = p, (j){t) = cj). The Metropolis 
algorithm generating initial configurations must respect the Gauss law. For this purpose the 
gauge-invariant Langevin equation of [4] was used in [5] to generate the initial {(/>, p}. The 
temperature dependence r(r) of these canonical measurements of the diffusion rate [5] is 
shown on Fig. I. At semiclassical temperatures (3 = > 7 the data agree well with the 
analytical predictions of [6], while at high temperatures /3 < 3 the rate exhibits the expected 
r^-law. In [7] it was found numerically advantageous to use Ag = - gauge in Cartezian 
coordinates during the Hamiltonian evolution. 

The gauge-invariant Langevin equation of [4] was introduced in the Abelian theory, where 
one can resolve the gauge constraints explicitly. To generalize it to the non- Abelian case we 
suggest to use the unitary gauge. In SU(2)-gauge theory with a scalar doublet the unitary 
gauge implies = \/ ip'^ip = p, and the Hamiltonian 

n = j-^E' + j-^H' + \p' + {Vpf + p'Aj + V{p) + ^ {D.Eif (5) 

can be expressed entirely in terms of canonically conjugated variables (/0,p) and {A1,Ef). 
The last term (Coulomb energy p'^A'^) takes into account the Gauss law DiEi = —2p^Ao to 
get rid of Ag. The Gauss law is automatically satisfied in, say. Metropolis samplings with the 
Hamiltonian (5). In this way the gauge-invariant Gibbs distribution in classical non- Abelian 
gauge theories can be performed by Monte-Carlo lattice simulations. Field configurations 
generated in the unitary gauge do not contain sphalerons. The latter appear during the 
canonical evolution according to the equations of motion, which are gauge-invariant. 
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Langevin and canonical times are two alternative definitions of physical time in the 
processes at T 7^ 0. To find a relation between them the sine- Gordon model was studied [8]: 



dx 



\ i^ f + \ { A^ f + 1 - cos{^) 



(6) 



The diffusing variable now is the zero mode (f) = 2^ / dx(f)(x). The model has many 
degenerate classical vacua and stable kink solutions interpolating between them. At T 7^ 
the Brownian motion of cj) has been observed. At semiclassical temperatures it is due to 
creation of kink-antikink pairs, so the thermal activation-type of the temperature dependence 
r ~ expl—Ex/T] has been seen in both Langevin and direct Gibbs measurements. 

The complete second-order Langevin equation has friction coefficient 7 as an input pa- 
rameter: 

9^ = -l(Px + W), (7) 

with the white noise normalized as < i]x(t) r]x(0) >= 2TjS(x)S(t). Eq. (7) interpolates 
between the Hamiltonian equation of motion for 7^0 and the first-order Langevin equation 
(no (^2;-term in (7)), valid at large 7 or t. The absolute value of the rate (the preexponential 
factor) depends on 7. For moderate to large friction the simulations of (7) show (Fig. 2): 



r(7) 



which would be a prediction of the first-order Langevin equation. That is 



because F is extracted from the correlator at large times for which the first-order Langevin 
equation is usually derived. 
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Figure 2: F ^ vs friction coefficient 7 at T = 2 in the sine- Gordon model [8]. 

In the limit 7^0 the rate remains nonzero and has been found [8] to coincide exactly 
with the corresponding independent direct Gibbs measurement according to (4). These data 
demonstrate a fundamental conclusion that canonical measurements can be obtained as the 
7 0-limit of the second-order Langevin simulations. 

In the high-temperature domain T > Ek at moderate to large friction the second- 
order Langevin simulations [8] show: F = 2T/j, which is obviously the prediction of the 
free first-order Fangevin equation. These simulations are insensitive to the dynamics of the 
system, because at high temperatures the noise-term i]x(t) entirely dominates over the regular 
Hamiltonian force. 

The high-temperature domain of canonical measurements was studied in [9]. It corre- 
sponds to the second-order Fangevin simulations with asymptotically small 7, so that the 
noise-term does not dominate over the regular force at any fixed T. This set up represents 
a large isolated system like the early Universe with no external dumping. In addition to the 
expected temperature behaviour F = a strong lattice spacing dependence of the diffu- 
sion rate was discovered (Fig. 3): F ~ . This is in contrast with the moderate to large 
friction simulations which do not exhibit any lattice-spacing dependence for a < 1 [8]. The 
rate, however, remains remarkably constant for the fixed ratio T/a even as one goes into the 
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Figure 3: r(r) for a = (1,0.5,0.25) at high T in the sine-Gordon model 



semiclassical temperature domain [9]. Presumably, damping with the friction 7 screens the 
high-frequency ultraviolet modes uj ~ > 7, which are essential beyond the semiclassical 
domain of temperatures. The lattice spacing dependence of the diffusion rate may be different 
between the sine-Gordon and Abelian Higgs models, because the dimensionality of the rate is 
different. The latter is fixed by the requirement that the diffusing variable is dimensionless, 
since it has integer dimensionless values in the classical vacua. Study of the lattice spacing de- 
pendence of the diffusion rate requires rather long simulations. As a decreases, the number of 
degrees of freedom grows, and it takes more and more time for the system to thermalize itself 
in the course of the canonical evolution, which always takes place before the Brownian motion 
sets in. This study is, however, important particularity in higher space-time dimensions, since 
it is a connection between the simulated classical system and the underlying quantum field 
theory at finite temperature [2], [10]. 
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